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ABSTRACT 

Pulse timing of the isolated pulsar PSR B1882-11 shows strong Fourier power 
at periods ~ 1000, 500 and 250 d, correlated with changes in the pulse profile 
(Stairs, Lyne & Shemar 2000). We study the extent to which these data can be 
explained by precession of the star's rigid crust coupled to the magnetic dipole 
torque. We find that the correlated changes in the pulse duration and spin period 
derivative can be explained as precession at a period of ~ 500 d with a wobble 
angle of ~ 3° if the star's dipole moment is nearly orthogonal to its symmetry 
axis. The dipole torque produces a harmonic at ~ 250 d. Comparison of the pre- 
dicted spin dynamics with the observed pulse durations requires the radio beam 
to have a non-standard "hour-glass" shape. We make predictions of variations 
in beam polarization and pulse profile with which to test this interpretation. 
The precession interpretation of PSR B1828-11 seriously challenges the current 
understanding of the liquid interior of the neutron star. In particular, if the in- 
ternal liquid is in a superfluid state, its rotational vortices cannot be significantly 
pinned to the crust. 

Subject headings: dense matter — magnetic fields — stars: magnetic fields - 
stars: neutron — pulsars: individual (PSR B1828-11) 
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Table 1: Spin Parameters of PSR B1828-ll a 



period p 


0.405 s 


V 


60.0 x 10~ 15 s s^ 1 


Spindown age t agc = p/2p 


0.11 Myr 



"Stairs, Lyne & Shemar 2000 

1. Introduction 

The monitoring of pulsar spin behavior through precision timing measurements pro- 
vides a means for studying the dynamics and internal structure of neutron stars. In this 
connection, detection of precession in pulsars is of particular interest. Relativistic (geode- 
tic) precession has been observed in the double neutron star binaries PSRs B1913+16 and 
B1534+12 (Weisberg, Romani & Taylor 1989; Arzoumanian 1995), and provides confirma- 
tion of General Relativistic orbital mechanics. Classical precession in isolated pulsars could 
probe the coupling of the stellar crust to the interior, but its detection is complicated by 
the presence of timing noise. Evidence for long-period precession is seen in the Crab pulsar 
(Lyne, Pritchard & Smith 1988), the Vela pulsar (Deshpande & McCulloch 1996) and PSR 
B1642-03 (Cordes 1993; Shabanova, Lyne & Urama 2001). The 35-d periodicity seen in the 
accreting system Her X-l (Tannanbaum et al. 1972) has been interpreted as precession by 
many authors {e.g., Brecher 1972; Trumper et al. 1986; Cadez, Galicic & Calvani, 1997; 
Shakura, Postnov & Prokhorov 1998). 

Recently, Stairs, Lyne & Shemar (2000) discovered the most compelling evidence to 
date for precession of an isolated neutron star. The period residuals of PSR B1828-11 
have an amplitude of ~ 1 ns and have strong Fourier power at periods of ~ 1000 d, 500 
d and 250 d. These Fourier components alone might suggest a Doppler effect through 
gravitational interaction with planets. However, correlations of the pulse shape and duration 
with variations in the period derivative p strongly suggest that the viewing angle with respect 
to the emission beam is varying, as one would expect if the star were precessing (see Fig. 
1). There is no obvious damping of the observed fluctuations, though the timing behavior 
is not entirely stable. See Table 1 for the spin parameters of PSR B1828-11. 

The purpose of this paper is to investigate if the timing behavior and pulse profile 
changes of PSR B1828-11 can be quantitatively understood in terms of stellar precession. 
We assume that the observed timing residuals represent true changes in the pulsar's rotation 
rate. An immediate question is which of the three observed Fourier components corresponds 



- 3- 



to the actual precession frequency. The residuals of the period p and period derivative p 
show the strongest Fourier power at ~ 500 d (see Fig. 1 and Stairs, Lyne & Shemar 2000 
Fig. 3); we tentatively assume that this Fourier component represents the precession period. 
The wobble angle 9 between the star's symmetry axis and its angular momentum can be 
estimated from the pulse profile changes. The difference in widths between the "wide" and 
"narrow" beam profiles is ~ 2° (see Fig. 2), implying a wobble angle of similar magnitude 
for plausible beam shapes. 

This precession interpretation could be complicated by magnetospheric processes. Pre- 
cession of the star could create changes in the magnetosphere that affect the geometry and 
location of the emission region. Changes in the height of the emission region would in- 
troduce light-travel time contributions to the pulse arrival times. Different regions of the 
magnetosphere might be viewed as the star precesses. We do not consider these effects in 
this paper. 

In the next section, we give an overview of precession. In §3 we study how the electro- 
magnetic torque, responsible for the star's spin down (at least in part), modifies precession. 
In §4, we relate this torque effect to the observed timing. In §5, we propose a model of 
the the observed pulse duration variations. In §6 we discuss some of our chief results. We 
conclude with discussion of observational tests of the precession interpretation. 

2. Overview 

Most treatments of pulsar precession to date have considered mainly free precession 
(e.g., Goldreich 1970; Pines & Shaham 1974; Shaham 1977; Suto & Iso 1985; Orford 1987; 
Bisnovatyi-Kogan, Mersov & Sheffer 1989; Glendenning 1990; Nelson, Finn & Wasserman 
1990; Bisnovatyi-Kogan & Kahabka 1993; Cadez, Galicic & Calvani 1997; Shakura, Postnov 
& Prokhorov 1998; Sedrakian, Wasserman & Cordes 1999; Jones & Andersson 2000). The 
observed variations in p and p in PSR B1828-11, however, are far too large to be explained 
by free precession for a wobble angle of 9 ~ 2°. In §3 we calculate the precession of a 
neutron star subject to a vacuum magnetic dipole braking torque and find a solution that is 
consistent with observations. 

Consider first the free precession of a biaxial neutron star; the geometry is as shown in 
Fig. 3. At any instant, the symmetry axis e$, angular momentum L and angular velocity 
u> span a plane. Let the oblateness be e = (/3 — h)/h, where /3 and l\ are the largest and 
smallest moments of inertia about the principal axes. For small oblateness e < 1, the angle 
between L and uj is 9 ~ e6 <C 9 (see, e.g., Jones & Andersson 2000). In the inertial frame, e 3 
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and u) rotate in a right-handed sense about L at frequency ~ uj{1 + e). Combined with this 
motion is a slow retrograde motion of the body about e 3 at precession frequency u p ~ euj. 

For a precessing pulsar, the pulse emission times are defined by the instants at which 
the beam center passes closest to the direction of the observer; the observed arrival times 
are simply offset from the emission times (apart from Doppler shifts from translational 
motion). If the star is freely precessing, the pulse arrival times vary sinusoidally as At ~ 
(po/27i)8 cot xsin(u;pt + where po is the spin period, x > $ is the inclination angle between 
the center of the pulsar beam and the body's symmetry axis and (3 is a phase (see Nelson, 
Finn & Wasserman 1990). This variation is a purely geometrical effect, as the magnitude of 
the angular velocity is unchanged by free precession. The corresponding changes in period 
derivative due to this geometrical effect are Ap ~ (2ix)~ 1 (pquo p ) 2 9 cot x sm(u p t + (3) (see eq. 
21 below). For a precession period of 500 d and x > 30°, the residuals in p are less than 0.02 
s s -1 , a factor of at least 10 smaller than observed. 

The electromagnetic torque on a neutron star modifies the precession and increases the 
magnitude of timing residuals for a given wobble angle (Jones 1988; Cordes 1993; Jones & 
Andersson 2000). In the body frame of a precessing object, the angular velocity vector u> 
moves about the body's symmetry axis e 3 at frequency oj p in a right-handed sense. For a 
magnetized neutron star with a magnetic moment misaligned with respect to the symmetry 
axis, the angle between the dipole moment m and uj changes by about 9 over a precession 
period (see eq. 22 below). The external dipole torque, which depends on this angle, thus 
varies over a precession period in addition to its slow decay over an evolutionary timescale. 
For small wobble angle 6, these torque changes produce changes in the star's spin rate u of 
(see eq. 13) 



For simplicity we assume that the pulsed radiation is emitted along the axis of the magnetic 
dipole, i.e., both the beam and the dipole are at an angle x with respect to the body's 
symmetry axis. The torque-induced timing variations are much larger than the geometrical 
effects associated free precession (Cordes 1993; Jones 1988; see eq. 21 below). In addition, 
and as pointed out by Jones & Andersson (2000), the dependence of the torque on the angle 
between the star's angular velocity and its dipole moment creates a strong harmonic of the 
precession frequency in the star's spin rate if the star is a nearly orthogonal rotator. In §5 
we apply this picture to PSR B1828-11 and find that x — 89° an d a wobble angle of 9 ~ 3° 
can account for the 500-d Fourier component with a harmonic at 250 d. Precession does not 
give the whole story, however, as it does not explain the Fourier power at 1000 d. 

As the star precesses we see different sweeps through the beam of the pulsar, allowing 
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us to partially map the beam morphology. As we will show, the observed variations in pulse 
duration imply that the beam has an "hour glass" shape, i.e., the pulse duration is longer 
for sweeps farthest from the beam center. 

3. Precession Under the Dipole Torque 

Precession of a neutron star will occur if the symmetry axis of the crust is not aligned 
with the angular momentum. The elastic crust has two contributions to its oblateness: a 
centrifugal deformation associated with the star's spin and a Coulomb deformation sustained 
by the rigidity of the lattice (Alpar & Pines 1985; Jones & Andersson 2000). The princi- 
pal axis of inertia of the centrifugal deformation follows the instantaneous spin axis of the 
crust, and so does not affect the precession. In §6 we discuss possible origins of Coulomb 
deformation. 

The period and amplitude of the crust's precession are affected by the coupling of 
the crust to the liquid interior. The precession creates time-dependent velocity differences 
between the crust and liquid that vary over the star's spin period. If the coupling time 
Tf between the crust and the liquid interior is much longer than the crust's spin period p, 
the precession will damp over ~ 2nTf/p precession periods (Bondi & Gold 1955; Sedrakian, 
Wasserman & Cordes 1999). All studies of damping of differential rotation between the crust 
and various parts of the liquid interior give r/ 3> p (e.g., Alpar & Sauls 1988; Epstein & Baym 
1992; Jones 1992; Abney, Epstein & Olinto 1996; Mendell 1998). To a good approximation, 
therefore, the liquid interior can be treated as decoupled from the solid (Bondi & Gold 1955; 
Sedrakian, Wasserman & Cordes 1999); in this regime the crust precesses almost as if the 
liquid interior were not there. 2 

In our calculations we treat the neutron star as consisting of a biaxial, rigid crust afloat 
on a liquid core. The solid portion of the star comprises less than 1% of the star's total 
moment of inertia and we regard the liquid of the inner crust as an extension of the core 
liquid. We neglect angular momentum exchange between the solid and the liquid. For 
simplicity we ignore the centrifugal bulge (which does not affect the precession) and treat 
the crust as perfectly rigid. For illustration we describe the spin dynamics of the crust with 



2 A precessing rigid shell containing a non-spherical normal fluid is inertially coupled to the fluid (Lamb 
1952). The fluid, which tends to a configuration symmetric about its rotation axis, exerts a reaction force on 
the solid. The relevance of this effect to an elastic neutron star containing superfluid neutrons and protons is 
unclear. Inertial coupling would not change the precession dynamics that we calculate in this paper, though 
it would change the dependence of the precession period on the crust oblateness. 
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the vacuum dipole torque of Davis & Goldstein (1970). While plasma effects could contribute 
to the spin down, a torque similar to the vacuum torque will always be present. In a frame 
corotating with the crust, the equations of motion are (Goldreich 1970; Melatos 2000) 

T T 2uj 2 , . 1 . 

I • — + u> x L = x m) x m + (w • m)(w x m), (2) 

where I is the moment of inertia tensor of the crust, uj is the angular velocity, L is the angular 
momentum, m is the magnetic dipole moment (centered on the star) and R is the radius of 
the star. The first contribution to the torque is due to the radiation far field and is responsible 
for spinning down the star. The second term represents the near-field radiation torque. The 
coefficient of this term depends on how the star's internal magnetization is treated (Melatos 
2000); its precise value does not affect our results. We will find that variations in the spin 
rate associated with precession are determined predominantly by the far-field torque (see. 
eq. 13). 

We define a Cartesian coordinate system (1,2,3) with the 3 axis along the symmetry 
axis e 3 of the crust. Without loss of generality, we take the magnetic dipole to lie in the 1-3 
plane, inclined with respect to the 3 axis by an angle x < 7r/2. At t — the spin rate is 
ujq = 2ir/p Q . The equations of motion are 

hui + u 2 u 3 (I 3 - h) = Ti = 7^-/3(^3 sin x-vi cos x) cos x 



2t 



3Xn 



H — 7^-/3(^1 sin x + <^3 cos x)o; 2 sin x (3) 
4r 



1 

hu 2 ~ U!UJ 3 (I 3 - h) = T 2 = -—I 3 u> 2 



2r 



3Xn 



+-^rh(wi sin X + ^3 cos x) (^3 sin x~wi cos x) (4) 

J 3 ^3 = T 3 = ^-I 3 sin x(ui cos x - ^3 sin x) 
It 

3x~ x 

— ^-/ 3 (^isinx + ^3Cosx)w2sinx, (5) 

where r = 3c 3 I 3 /4m 2 UQ is the characteristic time over which the far- field torque spins down 
the star, m = |m| and x = Ru /c « 1. For the oblate star, I 3 > I 2 = I x . Terms on 
the right-hand side of eqs. [3]- [5] proportional to Xq 1 arise from the near-field torque; the 
remaining terms are due to the far-field torque. The near-field torque is stronger than the 
far-field torque by a factor ~ x~ 1 , though its direction is always orthogonal to uj. 

Identification of the small parameters in the problem allows perturbative treatment of 
the torque. Because the torque acts over timescales that are much longer than the observed 
precession period, we expect the spin solution to resemble that for free precession. Hence, 
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the ratio uj p /uq, where uj p is the frequency of free precession, should be approximately e = 
(J 3 - h)/h. This ratio is ~ 10~ 8 for PSR B1828-11 if the precession period is ~ 500 d. 
Changes in the pulse duration (Fig. 2) through the precession period suggest a wobble angle 
of 9 ~ 2°. For these small values of e and 9, 9 ~ e9 <C 9, so # and a are nearly equal (Jones 
& Andersson 2000). In terms of the components of the angular velocity, a is approximately 
(cuf + ujI) 1 ! 2 / ujq <C 1. We therefore regard uj\ and uj 2 as perturbations about the secular spin 
down. The relative magnitudes of the small dimensionless parameters are given by 

ill 

1 > a > e > ea > x L > . (6) 

UJ G T u t 

To order eu 1 and eu 2 , eqs. [3] and [4] are 

o)i + eid 2 ^o = (7) 

u) 2 — ecJic^o = 0. (8) 

For a < 1, the solutions are 

ui = u a cos(cj p t + pi) (9) 

u>2 = u asm(u> p t + P) (10) 

l3 — h / ni N 

cj p = ecj = — 7 ^o, (11) 

h 

where P is a phase. The behavior of the component of the angular velocity orthogonal to 
the symmetry axis is essentially unaltered by the torque; the precession frequency is that of 
free precession and a ~ sin a = {uj\ + u^Y^/ujo =constant. In the body frame, both L and 
uj (which are always coplanar with e 3 ) precess slowly at frequency uo p uj about £3 and in 
the same sense. For the limit e C 1 that we are considering, a and 9 are nearly equal. 

For biaxial free precession, uj = is constant. Under the dipole torque, however, u 
undergoes changes that are important for the observed timing behavior. Combining eqs. 
[3]- [5] gives 

^ 2 C.T- e :W). (12) 



dt h V ^ h 

The uj ■ T term does not depend on the near-field torque. The near-field torque changes the 
spin rate only through the small final term, which we will neglect. The secular spindown in 
the absence of precession is given by uj ~ — sin 2 x/2r. Combining eq. [12] with eqs. [9]- 
[11], setting a = 9, retaining terms to order 9 2 {uiqt)~ 1 , and subtracting the secular spindown 
gives (dropping constant terms): 



Au _ 9 
uj 2t 



sin 2x cos(uj p t + P) + - sin 2 x cos(2uj p t + 2/5) 



(13) 
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Hence uj undergoes small variations of order 9(ujqt)' 1 and # 2 (cc;ot) _1 , driven by changes in 
the far-field dipole torque. These torque variations are due to changes in the angle between 
uj and m as w precesses through the body. The harmonic at twice the precession frequency 
arises from quadratic dependences of uj{T\ and uj-iT^ on the components of the angular 
velocity; the harmonic dominates the fundamental for 9 tan x > 4. For a small 9, the dipole 
moment must be nearly orthogonal to the rotation axis for the harmonic to be significant. 
The small oscillations of uj largely determine the observed pulse timing; they give period 
variations Ap ~ — (pl/2n) ACj. 



4. Timing 

The observed pulse timing at a given time in the precession cycle depends on the star's 
spin rate and the orientation of the angular velocity uj with respect to the observer. To 
relate the solution of the previous section to the observed pulse timing, we go to an inertial 
coordinate system x, y, z with the angular momentum along the z axis and the observer in 
the x — z plane with x > and z > 0. The azimuthal and polar angles of the magnetic 
dipole m, $ and 6, are given in this frame by 

^_ ^ m y (cos -0 sin + cos # cos sin -0) sinx — sin # cos cos x 
m x (cos -0 cos <fi — cos 9 sin sin ip) sin x + sin 6> sin cos x 

cos = — - = sin 9 sin ip sin x ~t~ cos 9 cos Xi (15) 
m 

where ip, the wobble angle 9 and are Euler angles following the definitions of Landau & 
Lifshitz (1976) and Goldstein (1980). 

The pulsar beam is not necessarily in the same direction as the dipole moment. Nev- 
ertheless, for simplicity we define a pulse as occurring when the azimuthal angle $ of the 
magnetic dipole equals the azimuth of the observer, i.e., when m y = and $ = 0. Hence $ 
is the observed pulse frequency. For small 9 < x an d #cotx 1, eqs. [14] and [15] give to 
second order in 9, 

9 2 

$ = + ?/> + 6> cot x sin ^ - — (1 + 2 cot 2 x) (2 cos 2 ^ - 1)^- (16) 
The angles ip and are given by 

rj} = tan" 1 ip- = tan' 1 — = ^ - uo p t - (3 (17) 

= ^-w + ^ + X- (18) 
cos 9 2 

(19) 
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The observed beam sweep rate is (dropping constant terms) 

$ = to - 9 uj p cot x cos(uj p t + 13)- 9 2 ujp (1 + 2 cot 2 x) cos 2 \uj p t + (3) (20) 

Hence, the sweep rate is the spin rate plus a modulation at the free precession frequency 
and its first harmonic. This modulation is a geometrical effect due to variations in the 
rate at which the beam sweeps past the observer as the symmetry axis precesses about the 
angular momentum; it is the only modulation that occurs if the precession is free. The sweep 
rate also changes as oo changes through secular spin-down and the variations Au; associated 
with torque variations as the star precesses. We henceforth focus on behavior of the period 
derivative since, as we will show, the torque changes evaluated in the previous section largely 
determine the observe spin behavior. Differentiating eq. [20] in time and combining with 
eqs. [13], [17] and [19], gives the period derivative residuals with respect to the secular spin 
down: 

q 

cot x cos(u p t + (3) + ~x cos(2cj p t + 2/3) 

- ^-(p LUp) 2 6 [cot x sin(cj p t + (3) + 9(1 + 2 cot 2 x) sin(2cu p t + 2(3)] . (21) 

2ir 

Here t age = rj sin 2 x is the characteristic spin-down age. This equation gives the complete 
expression for the star's residuals in p. The first set of square brackets gives the residuals 
due to torque variations from eq. [13]. The second set of square brackets gives the residuals 
due to geometrical effects (see Nelson, Finn & Wasserman 1990; Jones & Andersson 2000). 
For PSR B1828-11, the first-order torque effects always dominate the geometrical effects by 
a factor of ~ 250 if the precession period is ~ 500 d. For the harmonic at frequency 2lu p to 
play a role in the timing (terms proportional to 9 2 ), x must be nearly ir/2. The near-field 
torque, because it does not significantly change cu, is unimportant for the small-6 1 precession 
that appears to be taking place in PSR B1828-11. 

We model the data for the 2000-day span beginning at MJD 49,000 as the observations 
were most closely spaced during this period and these data show the observed periodicities 
most clearly. Fitting eq. [21] to the data, neglecting the small geometrical terms, gives the 
parameters 9, u p , (3 and x featured in Table 2. Comparison of the fit to the data is shown 
in Fig. 1. 

5. Pulse Duration Variations 

The duration of PSR B1828-ll's pulses is seen to change with the same periodic struc- 
ture as the spin residuals, suggesting different sweeps through the beam as the star precesses. 
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Table 2: Model Parameters 



e 3.2° 

X 89° 

2tt/u p 511 d 

f -0.3° 



To study this further, consider the behavior of the polar angle of the beam with respect 
to L (eq. 15). Expanding eq. [15] to first order in 9, using eq. [17], gives 

= x -6cos(uj p t + [3). (22) 

As the star precesses, G varies sinusoidally about x- The variations in the polar angle are 
in phase with the dominant contribution to the torque variations (the first term in eq. [21]). 

To quantify this further, suppose the constant angle between the observer and the 
angular momentum is 7 = £ + x ( see Fig. 4). We define the sweep angle AO as the 
difference in polar angle between the observer and the dipole at the time of the pulse: 
A© = £ + 9 cos(u p t + f3). As a simple model of the variations in pulse duration, we suppose 
that the pulse duration w is a function of A6 only. We assume that w has an extremum (it 
will turn out to be a minimum) for A6 = 0, corresponding to a viewing angle co-linear with 
the dipole axis. We take the expansion 

w = w + w 2 (AQ) 2 . (23) 

Consistent with Stairs, Lyne & Shemar (2000), we define a shape parameter S as 

w = max (w)S + min (w)(l — S), (24) 

where max (w) and min (w) are the maximum and minimum values of the beam duration 
in the precession cycle. Combining eqs. [22] and [23] gives 

S={ Vl«l+ e 7 25) 

l(l&) «*<°- 

This shape parameter depends on £, 9, u p , (3, and the sign of w 2 but not its magnitude. 

With 9, u p and (5 given by the fit to the timing data, the shape parameter is determined 
entirely by £ and the sign of w 2 - Interestingly, we find that w 2 must be positive, i.e., the 
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beam duration is larger for beam sweeps farthest from the dipole axis. Such a beam pattern 
is not standard, but is requiredby the precession interpretation of PSR B1828-11. The best- 
fit value of £ is —0.3°, as given in Table 2. Comparison of the fit to the data is shown in Fig. 
1. 



6. Discussion 

We find an acceptable fit to the data with x = 89° between the magnetic dipole and 
the star's symmetry axis, a precession period of 511 d and a wobble angle 6 = 3.2°. The 
predicted shape changes are in qualitative agreement with the observations if the angle 
between the observer and the angular momentum is 88.7° and the radio beam pattern has 
the "hourglass" shape sketched in Fig. 5. The beam shape that we have inferred differs 
form the usual picture of a roughly circular pulsar beam, but might be explainable in terms 
of a patchy beam of the type found by Han & Manchester (2000) in their statistical study 
of pulsar beam morphology. If the emission is beamed both parallel and anti-parallel to the 
dipole axis, an inner pulse would be expected for the nearly-orthogonal inclination that we 
have inferred. None is seen, so the emission pattern must be more complex. 

Our model does not account for 1000-d Fourier component seen in the data. However, 
significant aperiodic variations in the timing seen over timescales of years indicate that more 
than just precession is at work in PSR B1828-11. Slow changes in the emission region could 
be a contributing factor. Perturbations to the star's figure brought about by relaxation of 
the star's structure as it spins down could also play a role. Such starquakes would excite 
precession by perturbing the star's angular velocity with respect to the symmetry axis (Link, 
Franco & Epstein 1998). The perturbations would introduce irregular contributions to the 
magnitude of the wobble angle and the precessional phase (Jones 1988). 

Andersson & Jones (2000) have suggested that the 1000-d Fourier component represents 
the precession period, and that near orthogonality of the dipole moment with respect to the 
symmetry axis gives a harmonic at 500 d in the spin residuals. This scenario is unlikely for 
several reasons: 1) near perfect orthogonality (x~89.9°) is required to keep the first-order 
torque term in eq. [21] from giving too large a contribution to p, 2) the strong Fourier 
component at 250 d is not accounted for, and, 3) the shape parameter, which would now be 
periodic at 1000 and 500 d, would bear little resemblance to the data. 

The stellar oblateness inferred from the precession period of 511 d is e = (I3 — h)/h = 
9 x 10~ 9 . For comparison, we estimate the maximum Coulomb oblateness. The maximum 
mass M m of a "mountain" that the rigid neutron star crust can support is given roughly 
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by M m g ~ B6 C / R where g is the gravitational acceleration at the stellar surface, B is the 
bulk modulus of the crust (~ 10 48 erg) and 9 C is the critical strain angle at which the crust 
fractures. The maximum oblateness sustainable by crust rigidity is thus (see also Ushomirsky, 
Cutler & Bildsten 2000) 

e = — ^wr l{] 10=*' (26) 

where M c is the mass of the crust. A mountain could form as the result of a starquake that 
perturbs the star's principal axis of inertia away from its angular momentum axis (Link, 
Franco & Epstein 1998). However, even if the crust is relaxed (i.e., unstrained) before it is set 
into precession, rigidity prevents a portion of the star's bulge from following the instantaneous 
rotation axis. The effective oblateness due to rigidity is in this case e r ~ 10~ 5 en, where is 
the centrifugal oblateness (Alpar & Pines 1985; see also Munk & MacDonald 1960). For PSR 
B1828-ll's spin period e r ~ 2 x 10 _11 (p/0.4 s) -2 , giving a precession period of 700 yr. We 
conclude that PSR B1828-ll's crust is strained to the extent that the Coulomb oblateness 
is ~ 100 times larger than e r . A significantly strained crust can precess at a much higher 
frequency than the unstrained crust considered in some previous studies (e.g., Goldreich 
1970 and Melatos 2000). 

The precession of PSR B1828-11 seriously challenges the notion that superfluid vortices 
of the crust pin to nuclei (e.g., Anderson & Itoh 1975). Pinning would drastically alter 
the way in which the star precesses. Shaham (1977) showed that very effective pinning 
changes the precession period to (Ii/I p )p, where I p is the portion of the star's fluid moment 
of inertia that is pinned and Ii is the moment of inertia of the part of the star that is 
precessing. Sedrakian, Wasserman & Cordes (1999) have shown that Shaham's conclusion 
is essentially unaltered if the pinning is imperfect. Pulsar glitches, which might arise from 
variable coupling between the superfluid and the crust, indicate that / p // s tar > 1.4% in stars 
that frequently glitch (Link, Epstein & Lattimer 1999); / sta r is the total moment of inertia 
of the star. This degree of pinning in PSR B1828-11 would give a precession period of <C 40 
s, far shorter than the observed of 500 d precession period. It may be that vortex pinning is 
inhibited in stars that undergo precession with amplitudes as large as those in PSR B1828-11. 



7. Observational Tests 

Our proposed model for the fluctuations in the timing residuals of PSR B1828-11 implies 
that the pulsed emission should have several specific attributes that may be observationally 
tested. Fig. 5 illustrates these features that are characteristic of our model. During the 
511 day precession cycle the observer sees both the broad upper and lower parts of the 
beam. One important test of our model would be evidence of these two different parts of 



-13- 



the emission region. 

We modeled the observed pulse durations with the parabolic form of eq. (23), giving the 
overall beam shape illustrated by the dashed lines in Fig. 5. The actual shape of the emission 
must be more complex, however, as we suggest with the irregular brightness contours also 
shown in the figure. Additionally, the linear polarization of the emission region is unlikely 
to be the same in the upper and lower portions of the emission region, as indicated by the 
double- headed arrows in the figure. 

To indicate the possible signatures one may observe, we have marked certain viewing 
angles (A-C) in Fig. 5 and the corresponding observing times in Fig. 1. The dotted 
horizontal line in Fig. 5 is the shape equator (perpendicular to the largest principal moment). 
The horizontal dotted lines A and C give the range of the viewing angles seen from Earth. 
Line B represents the viewing angle at which the pulse profile is the narrowest. The pulse 
profile is wider at viewing angle A and widest at viewing angle C. 

While the shape parameter is the same at phases above and below phase B, the precise 
pulse profiles and radio polarization should be noticeably different at these phases. The 
polarization angle will vary across the emission region so that the polarization direction will 
change continuously from phase A to C and then reverse going from C to A. Detection of 
these changes in pulse shape and polarization would indicate that the radio beam has the 
hour-glass shape we have inferred, and would support the interpretation that PSR B1828-11 
is a precessing neutron star. 

We thank A. G. Lyne, D. I. Jones, R. Bandiera, A. Melatos, J. Cordes, D. Chernoff and 
the other participants of the ITP workshop on Spin and Magnetism in Young Neutron Stars 
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auspices of the U.S. Department of Energy, and was supported by NASA ATP grant NAG 
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Fig. I. — Timing and beam shape data for PSR B1828-11 (from Stairs, Lyne & Shemar 
2000). Only the most-densely sampled portion of the 13- year data span is shown. The top 
panel gives the period residuals with respect to the star's secular spin down. The middle 
panel is the time derivative. The bottom panel shows the shape parameter of Stairs, Lyne 
& Shemar (2000); it is S — An/(An + Aw), where An and Aw are the fitted heights of the 
narrower and wider profiles, so that S ~ for wide pulses and S ~ 1 for narrow ones. The 
data points give average values of S, obtained by averaging S over multiple bins. The solid 
curves for Ap and Ap are the fit described in the text. The jagged line in the lower panel is 
a guide to the eye. 
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Fig. 2. — Templates for the wide and narrow beam profiles of PSR B1828-11 (courtesy I. 
Stairs). 
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Fig. 3. — Geometry of free precession of a biaxial neutron star in the inertial frame. The 
body's symmetry axis is denoted by £3, the angular momentum by L and the angular velocity 
by uj; the three vectors always span a plane as shown. The angles 9, 9 and a are constant, 
with 9 ~ e9 <C 9 for small oblateness e. The vectors e 3 and u> rotate about L at nearly the 
spin frequency u>. A dipole moment m fixed in the body and taking an angle x with respect 
to e 3 , rotates in a retrograde sense about £3 at frequency~ ecu. 
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Fig. 4. — Observing geometry. At the instant the dipole moment m is in the plane containing 
the angular momentum L and the observer, the angles are defined as shown. 
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Fig. 5. — The beam pattern for the radio emission from PSR B1828-11. The dashed parabo- 
las represent the analytic fit to the beam width, and the contour curves are one possible 
realization of the actual beam shape. As the star precesses, observer's line of sight varies by 
±6* about the star's shape equator, the dotted horizontal line. The solid lines A and C give 
the range of the viewing angles. Line B represents the viewing angle at which the beam is 
the narrowest. The double-headed arrows represent plausible polarization directions. 



